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Abstract. The field-angular dependence and anisotropy of the critical current 
density in iron-based superconductors is evaluated using a phenomenological approach 
featuring distinct anisotropy factors for the penetration depth and the coherence 
length. Both the weak collective pinning limit, and the strong pinning limit relevant for 
iron-based superconductors at low magnetic fields are considered. It is found that in 
the more anisotropic materials, such as SmFeAsO and NdFeAsO, the field-angular 
dependence is completely dominated by the coherence-length (upper-critical field) 
anisotropy, thereby explaining recent results on the critical current in these materials. 
In less anisotropic superconductors, strong pinning can lead to an apparent inversion 
of the anisotropy. Finally, it is shown that, under all circumstances, the ratio of c-axis 
and afe-plane critical current densities for magnetic field along the a&-plane directly 
yields the coherence length anisotropy factor s^. 
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1. Introduction 

The multi-band nature of superconductivity in iron-based superconductors has now 
been weU estabhshed [H While direct measurement of the magnitude of the 
superconducting gap on the different Fermi-surface sheets has been and is the most 
convincing method in this respect, it is not the easiest. Therefore, many researchers 
have resorted to the measurement of the anisotropy of the superconducting parameters 
such as the coherence length and the penetration depth 0,^[^I3[71IHI^ to prove the 
multi-band nature of superconductivity of these layered materials. This manifests itself 
through the inequality of the anisotropy factors = Cc/Cab and £x = Xah/^o which are 
the respective ratios of the coherence lengths and ^ah parallel and perpendicular to the 
anisotropy (c)-axis, and of the penetration depths A^^ and Ac for supercurrents running 
in the a6-plane and along the c-axis respectively. The anisotropy of the coherence length 
is usually obtained from that of the upper critical fields and for field oriented 
parallel to, and perpendicularly to the c-axis [31 (U EJ [6[ [8[ [ini E]. The anisotropy of the 
penetration depth is most often extracted from the anisotropy of the lower critical fields 
B^^ = ($o/47rA^^) In and = ($o/47rAa6Ac) In /^c parallel and perpendicular to the 
c— axis O El [12], or from direct measurements using microwave techniques [131 El ES] 
(here K.ab = Kb/^ab while = V^abK/^ab^c) 116^. The inequality of 6:^ and Sx was 
first observed [TTl UHl [El [20| , and theoretically described [21j in the archetypical two- 
band superconductor MgB2. It was shown that in the case of weakly coupled two- 
band superconductors such as MgB2, Ginzburg-Landau theory, which, for single-band 
superconductors, introduces the anisotropy ratio s = = Sx diS the square-root of the 
ratio of the effective masses rUab and rric perpendicular and parallel to the c-axis, is 
inappropriate [22j. 

More recently, critical current density measurements have been used to address 
the anisotropy of superconducting parameters in the iron-based superconductors. In 
particular, Kidzun et al [23\ and Hanisch et al [24] applied the anisotropic Ginzburg- 
Landau scaling procedure of Ref. f32], developed for single band superconductors, to the 
critical current density of epitaxial LaFeAsO and Ba(Fei_a^Coa;)2As2 films. They found 
that the extracted temperature dependence of the anisotropy parameter, ^(T), matches 
that of the anisotropy parameter determined from scaling of the upper critical field [2l] . 
This appears to validate the use of critical current density measurements for obtaining 
reliable values of the anisotropy parameters. This approach was taken further by Moll 
et al [25], who specifically fashioned micrometer-sized SmFeAsO bridges to measure 
the three independent critical current densities, for current running in the a6-plane 
and field oriented along the c-axis, and j^^ and j^^ for current in the a6-plane and along 
the c-axis, respectively, and field in the a6-plane (see the definitions in Fig. [T]). The 
three critical currents were also measured independently in single crystalline LiFeAs by 
using a Hall-probe array magnetometry technique [26]. These measurements allowed 
the direct extraction of the coherence length anisotropy. 

However, the interpretation of the critical current anisotropy is, a priori , not so 
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Figure 1. Sketch defining ttie three independent critical current densities in an 
anisotropic superconductor, (a) j^^ is the (most commonly measured) critical current 
density for current in the a6-plane, field along the c-axis. (b) the a6-plane critical 
current density j^^ for field in the a6-plane. This configuration involves the hard 
motion of vortex lines along the c-axis. (c) the c-axis critical current density j^^ for 
field parallel to the a6-plane. This involves the easy motion of vortices along the 
a6-plane. 



straightforward. In principle, the critical current depends in a non-trivial manner on 
the specific nature of the pinning potential, that is, on the kind of defects and impurities 
that are present in the superconducting material. It also depends non-trivially on the 
details of vortex elasticity. A supplementary complication comes from the fact that 
vortex pinning in the iron-based superconductors is, certainly at the low magnetic fields, 
dominated by the so-called strong pinning ^271 [28] by nm-sized impurities ^29i [30l [3l] • 
In this regime, the anisotropy of the critical current may be influenced by the size and 
shape of the strongly pinning defects. Therefore, it is important to understand why 
the measurements of Refs. [23l El], as well as those of Konczykowski et al. [26], yield 
the coherence length anisotropy with such precision, even though the authors [231 [24] 
applied a single-band approach to analyze their results. 

In the following, we explore the effect of different anisotropics on the critical current 
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in the strong pinning regime in single-band and multi-band superconductors. Our goal 
will be limited to the identification of the main origin of the anisotropy of the critical 
current. This justifies a purely phenomenological approach, in which we will suppose 
that, in a single-band superconductor, the angular dependence of the properties of the 
superconducting state enter through the angular dependence of the penetration depth, 
A(i9) = \ah/^^ for supercurrents running at an angle i9 with respect to the a6-plane, 
while the angle-dependent coherence length ^{d) = ^ab^^- Here 1} is the angle between 
the orientation of interest and the a6-plane. The scaling function 

£^ = V (3Qg2 ^ _^ gij^2 ^ ^ 2 (1) 

can be found in Refs. [32l |33]. For multi-band superconductors, we will assume that 
the angular dependence of vortex-related quantities in the mixed state can be described 
by two different anisotropy ratios ^a('^) and s^{i}) for the penetration depth and the 
coherence length, and that the angular dependence of these ratios is the same as Eq. ([T]). 
In all cases, s^^ will be assumed to be larger than the distance between the structural 
layers, so that the material can be described as "continuously anisotropic". Josephson 
tunneling between layers is beyond our consideration. In section 2^ we first derive 
anew the results for the critical current density of Ref. ^33j, for the case of weak 
collective pinning of individual vortices, keeping track of the anisotropics introduced 



by the coherence length and the penetration depth, respectively. In section [2731 we then 
apply the same procedure to the formulae of Ref. [28j for strong pinning by large point 
defects, which is more appropriate in iron-based superconductors at low fields. Here, 
the effect of a non-spherical shape of the pinning defects is also taken into account. 

The procedure has the merit of demonstrating that the field-orientation-dependence 
of the a6-plane critical current density jabi"^) is, in all cases, determined by the 
competition of two factors. The coherence length anisotropy is responsible for the 
increase of the jabi'^) as the magnetic field is aligned towards the a6-plane, mainly 
because of the larger elementary pinning force experience by the vortex lines as these 
are turned towards the a6-plane. For large defects, this effect is counteracted by the 
penetration-depth anisotropy, which result in a larger line energy and therefore larger 
pinning energies when the vortices are aligned parallel to the c-axis. In most situations, 
the latter effect is much weaker than the first, so that the anisotropy extracted from 
field- angular dependent measurements of the critical current yields 6:^. 

The latter is also true for the anisotropy of the critical current density itself. On 
very general grounds, it is shown that the ratio jf^/j^ is, in all cases, equal to the 
coherence length anisotropy ratio e^. This fact was recently exploited to extract the field- 
dependence of 6:^ in single crystalline LiFeAs t26j. We note that the results below are 
relevant not only for iron-based superconductors, but also for high Tc superconducting 
cuprate thin films, tapes, and composite conductors. 
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2. Critical current anisotropy in uniaxial anisotropic superconductors 
2.1. Weak collective pinning in single-band superconductors 

The anisotropy of the critical current density may in principle, depend on the underlying 
pinning mechanism in the material. A flagrant example is the presence of anisotropic 
pinning centers such as heavy-ion irradiation-induced amorphous columnar defects 
[Ml ESI [36] , which introduce a preferential orientation for the vortex lines supplementary 
to that imposed by the anisotropic material. However, in the case of weak collective 
pinning by dense atomic-scale point defects, the anisotropy of the critical current density 
ought not to depend on the details of the defect landscape, but only on the material 
properties. It can thus be straightforwardly obtained using the scaling formalism 
introduced by Blatter et al [32l [33] . 

In the so-called single-vortex pinning regime, in which the critical current is not 
limited by intervortex interactions, the critical current density jab in the (ab) plane, 
perpendicular to the symmetry (c)-axis, is expected to be fleld- angle independent. 



Here, fp is the elementary pinning force of a single point defect for a vortex directed 
along the c-axis, Ud ^ is the point defect density, and the factors (^^^ and ^abCc derive 
from the statistical averaging of the pinning forces over the vortex core. The brackets 
(. . .) denote averaging over the vortex core. For fleld aligned along a6, one is faced with 
a larger pinning force for "out-of-plane" vortex motion along c, since the pinning 
range [38] rf ^ = ^Cab is shorter in this direction. However, this is compensated by 
the larger longitudinal correlation length Lf^ = Lc/s for fleld oriented along the a6-plane 
as compared to Lc for fleld || c ^3]. In this speciflc situation, Lf^ does not depend on the 
direction of the vortex displacement u in the pin potential. The anisotropic vortex line 
tensions, = e\je^ and = e\je for out-of-plane and in-plane vortex displacements 
respectively, compensate for the anisotropics of the elementary pinning force and the 
pinning range, i.e. 

( ^ v^^ 

L'c = (u ± ab) (3) 




(u II ah) (4) 

= Ule. (5) 

Here e\ ^ e^e^ and 6:0 = $o/47r/xoA^^ are the vortex line tension and the vortex line 
energy for magnetic induction B || c, and $0 is the flux quantum. Hence, the critical 
current density for fleld along the c-axis is equal to the critical current density 
(with B _L j) for fleld along the a6-plane. 
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The c-axis critical current density for B || ab on the other hand involves "in- 
plane" vortex motion, and is reduced by the anisotropy factor ^, such that j^^ = ej^^. 
These relations do not hold in the so-called "bundle regime", in which vortex interactions 
reduce the magnitude of the critical current density. Then, jab{B£^) follows the well- 
known field-angular scaling relation in = (^^ cos^ i9 + sin ^y/^ as detailed in Ref. [33j. 
Here i9 is the angle between the direction of the magnetic induction and the a6-plane. 



2.2. Weak collective pinning in multi-band materials 

We now turn our attention to the situation in the multi-band iron-based 
superconductors. In what follows, the multi-band nature of superconductivity will be 
taken into account phenomenogically, by simply introducing different anisotropics 6:^ and 
£x of ^ and A, respectively. Then we can rederive the results or Ref. [33], but keeping 
track of these different factors. We shall restrict ourselves to the single vortex regime of 
pinning, in which the critical current density is constant as function of magnetic field. 
In iron-based superconductors, this is thought to be relevant in a range of intermediate 
field strengths, sufficiently high for the critical current not be dominated by the strong 



pinning contribution due to nanometer-sized point-like pins (section 2.3), but below the 
onset of the so-called "second magnetization peak" regime, which is not described by the 
collective pinning theory [29l [30] . Typically, this corresponds to a field range of several 
tenths of Tesla, to a few Tesla at most. 



2.2.1. Elementary pinning force The elementary pinning force fp can, within factors, 
be written as the product of the condensation energy = £:o/27r^^^, the 

defect volume ^ttD^^ and the inverse range of the pinning potential, rj^. In the 
case of pinning by mean- free path variations ^Hl [301 1331 [37], a scattering factor 
proportional to g{pD){^o/ Dy) should be introduced, with g{pD) the Gor'kov function, 
Pd = hvF/27TTcl ^ ^o/l the disorder parameter, vf the Fermi velocity, / the mean 
free path, and = 1-35^(0) the temperature-independent BCS coherence length. 
Disregarding the possibility of anisotropic scattering, the elementary pinning force for 
magnetic fields not aligned with the c-axis is modified because of the anisotropy of the 
vortex line energy 6:0, which involves the kinetic energy of current components fiowing 
both parallel to c and ab as well as the anisotropy of the vortex core radius, and because 
of the anisotropic range of the pinning potential. When B is aligned along the a6-plane, 
the first effect introduces a factor Sx and the second a factor e^^. As for the pinning 
range, this also depends on the orientation of the supercurrent. For j || ab, j _L B, 
vortex motion is parallel to the symmetry axis and, in the low-field limit {B < 0.2S),'2), 
rj- = <^C5 while for j || c, = (^^6- Therefore, 

41=^/, (B||a5;j||c;|B|<0.25ll2) (6) 

= (B II ab;i II ab;i ± B; |B| < 0.2BI,). (7) 
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For higher fields, the pinning range is given by the intervortex distance: rj- = ao 
while rj = ao/ej^^ [with ao = (2$o/\/35)i/2 ^j^g vortex spacing for B 11 c 1 
Therefore 

41 = ^/,, (B||a6;j||c;|B|>0.2i?;2) (8) 



while 



fp = (B II ab;i \\ ab;j ± B; |B| > 0.2BI,). (9) 



'P 3/2' 



2.2.2. Pinning correlation length The isotropy of the pinning correlation length for 



B II ab in Eqs. (3|4) is nominally lifted, because the anisotropy of the vortex line tension, 
£i = ei/e\ and S\je\^ bears mainly to that of A, while the range r/ is determined 
by the value of the coherence length in the direction of the vortex displacement. As a 
result, the pinning length iX = {e'^'^ / e]l'^)Lc of a vortex confined to the a6-plane will 
no longer be equal to = [e^^'^ / e^^'^)Lc = {e^/ exY^'^L\ for a vortex confined to a plane 
containing the c-axis. 

However, in a three-dimensional superconductor, both in-plane and out-of-plane 
random displacements are allowed, and both determine the metastable vortex position 



and the critical current density. Taking into account Eqs. (6|7) and the difference of 



the line tensions for in-plane and out-of-plane vortex displacements, and minimizing the 
sum of elastic and pinning energy e\{^ah/Llf + e^{^c/Llf - {nd{fi^) / llf'^^^^CJ'^ - 
K(//')/4)^/'ec/'eir, we obtain 

22/3^1/3 

Ll = L/3 ^c. (10) 

If the coherence length anisotropy dominates the penetration depth anisotropy, the 
vortex excursions will be mainly confined to the a6-plane, the large pinning force in the 
c-direction prohibiting out-of-plane excursions. 

2.2.3. Critical current density As in Eq. ([2]), the critical current density is determined 
from the equality of the Lorentz force on a correlated pinned vortex segment and the 
total pinning force on that segment, 

%jcL,=(^'^y\c^^F, (11) 

where ^ is the square-root of the vortex core area. For the two different directions of 
motion of vortices parallel to the a6-plane, we find 

fc' = "-^^^^$^Jsv (B||a6;j||c) (12) 
fab = 5?3 (B||a6;j||a6;j±B). (13) 
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Note that for Sx ^ s^^ j^^ ejsv and j^l jsv-> in agreement with the results 
outhned in Section |2.1[ More importantly, the critical current density is independent 
of the range r/ of the pin potential. Since the core area and the pinning length are 
the same, irrespective of the direction in which the vortices are driven, the difference 
between these expressions is entirely due to the different elementary pinning forces for 
in-plane and out-of-plane vortex motion. As a result, the anisotropy 

■ah 

= (14) 

Jab 

is uniquely determined by the anisotropy of the vortex core {i.e. of the coherence length). 
Therefore, the measurement of the j^^/j^^-anisotropy directly yields the coherence 
length {Bc2-) anisotropy. 

The above conclusion, a consequence of the fact that the critical current density in 
the single vortex regime can be written as jc = Fp/^oLc^ can be extended to all field 
and temperature ranges. In particular, in the bundle pinning regime [33j, Eq. (11) is 
replaced by 

jcBV, = {n,{f,)Vcf' (15) 

where the correlation volume Vc = R^Lc is the product of the square of the transverse 
correlation length Rc and the longitudinal correlation length Lc. Since, for a given 
magnetic field orientation, the correlation volume is the same regardless of the direction 
in which the vortices are being driven, the anisotropy of the critical currents for different 
driving directions is always given by the anisotropy of the elementary pinning force. 



2.3. Strong pinning in anisotropic superconductors 

We now investigate the case of so-called strong pinning [271 [281 [39] , ^-^-^ the situation in 
which a type-II superconductor contains point-like defects or heterogeneities with typical 
dimensions 6^, parallel and perpendicular to the field direction larger the coherence 
length, and density much smaller than Then, the (first moment of the) pinning 
force from different defects does not average out as in the weak pinning case. The critical 
current density is now determined by the direct sum of the forces exerted by each pin 
a vortex line is able to take advantage of [27]. This situation is encountered in cuprate 
superconductor thin films [281 SO] and coated conductors [41] , in which second-phase 
nanoparticles are responsible for pinning, as well as in most iron-based superconductors 
at low fields [291 |30l [3l] • I^i the latter materials, the presence of nm-scale heterogeneity 
of the superconducting properties is thought to play the role of the strong pins, thereby 
determining both the highly inhomogeneous vortex distributions, as well as the low- field 
critical current density [31]. 

At low fields, the number of pins a vortex line can benefit from is solely determined 
by the vortex line tension. The mean characteristic length C = 

vortex between two effective pins is determined by the balance of the tilt deformation 
energy ei{u/CY of the line, and the pinning energy gain per pin. Up. At magnetic 
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inductions exceeding 5* = 7r$o^i (Up/so)^ the intervortex repulsion rather than the hne 
tension hmits the number of defects any given flux hne can take advantage of. The 
pinning length is then given hy C = riiU^^ where the trapping area = Up/ (c66C44)'^/^ao 
is determined by the balance between the pinning energy gain and the energy of the 
lattice deformation. For B || c, the shear modulus is given by Cqq ^ So/Aal^ and the 
non-local tilt modulus is evaluated as C44 ^ e'^So/al. For fleld along the c-axis, the 
critical current density is [28j 

where the maximum force that a single strong pin can exert on a vortex line (for B || c) 

is the product of the condensation energy over the range of the pinning potential, the 
core volume occupied by the defect, and a logarithmic factor taking into account the 
modiflcation of the supercurrent distribution around the defect [28[ . Here we have 
assumed that the pinning defects are insulating, and of ellipsoidal shape with half- 



axes bx \\ ab and 6^ || c. Eq. (18) interpolates between the logarithmic dependence on 
dimensions expected for large defects with bx^b^ ^ ^ab^ and the simple factor {bx/^abY 
proportional to the defect cross-sectional area that one has for small defects. 

2.3.1. Elementary pinning force In order to evaluate the anisotropy in the limit of 
strong pinning, both with respect to the fleld orientation and the driving direction, we 
again start by evaluating the change of the elementary force / with fleld orientation. 
Contrary to the case of weak pinning, / is reduced by a factor Sx for B || a6, due to 
the decrease of the vortex line energy. A supplementary complication enters in that 
the smaller core cross-section for fleld not aligned with c, as well as the geometrical 
anisotropy of the defects, affect the factor under the logarithm. For B || a6 and j || c, 
i.e. in-plane vortex motion, the pinning range r/ = ^ab is the same as for vortices 
parallel to the c-axis, while for the motion along c of vortices parallel to the a6, it is the 
smaller c-axis coherence length that enters. One thus has 

/" = ^/ (B II ab;} II c) (19) 

= f^^l /; (B II ab;3 \\ ab;} ± B) (20) 



for arbitrary orientation of the field with respect to the c-axis, 

/' ^ ll^/ (J II c) (21) 

f'-{7^))f- (JII-'^^J^B) (22) 
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where we have introduced the geometrical anisotropy 

1 In (1 + bl/2ej 



e,{d) = -= 

^Jsm' + I cos2 i3ln^l + b^^b'i sin' d + hi cos' ^/2Qf^e^{^) 

' (23) 

for model ellipsoidal defects. For spherical defects of size ^ ^ab^ this reduces to a simple 
factor s^{'d). The pinning energy for B || ab is 11^^ = [£x/£i){0)]Up^ while for arbitrary 
orientation it is Up{i}) = [£x{i}) /£i){i})]Up. 

2.3.2. Single vortex limit of strong pinning We first examine the limit of low fields. 
The pinning length C for arbitrary field orientation takes on different values depending 
on the direction of the deformation: it is larger for out-of-plane deformations, due to 
the larger line tension £^ ^ £i/£\{{}) in this direction ^33j. At equilibrium, the most 
favourable pinned configuration will be determined by 

- A (u II a,) (24) 

since this is smaller than 

— I p-^/^f??") — 

by a factor £x{'d) = ^J £\ cos^ d + sin^ d. Thus, the vortex adapts to the pin potential by 
deformations parallel to the ab plane rather than perpendicularly to the field direction. 
Note that in all cases, except that of extremely oblate defects (6:5 < 6:^), the typical 
length of a pinned segment || ab exceeds that for B || c, and thus, each vortex ought to 
be, on average, pinned by a smaller number of defects. If one assumes that the vortices 
can seek out the most favorable pinned configuration at each stage of their motion, the 
a6-plane critical current density simply follows as 

= ( ) = J'!:^\,, Js; (j II ab;i±B;B< B;) (26) 

valid for B < 5|* (see below). If the relevant pinning centers are relatively small (radius 
of the order of ^ab) and isotropic, the critical current density scales as 6:^(i9)^^^^^(i9). In 
the case of single-band superconductors, where £x{'d) = £^{'&) = 6:^, this reduces to the 
scaling law jabi"^) oc 6:^"^^^. 

For B II a6, the critical current densities parallel and perpendicular to the c-axis 
depend on the anisotropics of the penetration depth and the coherence length in a 



different combination than in the weak pinning limit of section (2.2): 



jf = 1 






Jab — 1 




1 _ 4(0) 







(j||c;S<S|p (27) 

j,; {3\\ab;}±B;B <Bl) 

(28) 
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Figure 2. (Color online) Angular dependence of the critical current density in 
the single vortex limit of strong pinning, for anisotropics representative of (a,b) the 
(Ba,K)Fe2As2 system and (c,d) the NdFeAsO and SmFeAsO systems, containing 
spherical insulating defects. The different values of are introduced to mimic the 
temperature dependence of the anisotropy. The graphs (a,c) depict the expected 
behaviour for small defect sizes, bx = bz = ^, while panels (b,d) are for large defects 
with bx = bz = 10^,. i9 is the angle with respect to the a6-plane. 



However, the generality of expression (16) together with a C that does not depend 
on driving direction signifies that their ratio only depends on the coherence length 
anisotropy. In particular, for B || a6, the ratio of critical currents for in-plane and 
out-of-plane motion of vortices 

jab 

= (29) 

Jab 

again directly yields the coherence length anisotropy. 

The case in which vortices cannot equilibrate during driven motion was explored 
in Ref. [42j. Then, vortices can only take advantage of those pins that are located 
within a distance Ut of the plane of motion. The "dynamic" critical current density 
jtr = f /^oLfoT driven motion depends on the vortex length L trapped between adjacent 
pins. Analysis yields, for the critical current density for field along c, \42\ 



/TT \ 5/9 

J^r = nff[fJ C'"^ (30) 



Anisotropy of strong pinning in multi-band superconductors 



12 




Figure 3. (Color online) Angular dependence of the critical current density in 
the single vortex limit of strong pinning, for anisotropics representative of (a) the 
(Ba,K)Fe2As2 system and (b) the NdFeAsO and SmFeAsO systems, for different radii 
bx = bz of spherical insulating defects, is the angle with respect to the a6-plane. 



The a6-plane critical current density for arbitrary field direction is obtained as using 
the parameters corresponding to (hard) motion, perpendicular to the field direction. 
Inserting relevant parameters, we obtain 



.29/9 



5/9 



1 



Jtr 



(31) 

(32) 



In the case of small spherical defects and similar anisotropics for the coherence length 
and the penetration depth, jab,tr{'&) ^ ^"^J^ jtr- 
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For field along the a6-plane, the critical current densities along c and ab are 
determined by hard and easy motion respectively, so that 

^19/9 

^c,tr = 14/9 2/9 -^'^^ ^^^^ 



^29/9 

^ab,tr = 14/9 11/9^^'' ' ^^^^ 

Unlike the above results, the ratio of the "dynamic critical currents" jc\r/jabtr ^ ^c/^a°^^ 
does not uniquely depend on the coherence length anisotropy. Surprisingly, the derived 
dynamic critical current densities for vortex motion along the hard- and the easy 
direction are nearly equal in single band materials, because the anisotropy of the 
elementary pinning force cancels that of the vortex line tension. 

2.3.3. Three-dimensional strong pinning Turning to higher fields, our task is 
to determine the trapping radii Uq and Uq for in-plane and out-of-plane vortex 
displacements, respectively. These are determined by the shear- and tilt moduli 
for vortex lattice deformations within the ab plane, Cqq = £\{i})cqq and C44 ~ £44, 
and out-of-plane vortex deformations, Cqq = CQQ/ex{d) and ^ (^a(^)/^a)^ £44 , 
respectively [43l HI]. Here 9 = | — i9 is the angle with respect to the c-axis, and 

ex{6) = ^/ e\ sin^ 6 + cos^ 6 |33]. Taking into account the distortion of the vortex lattice 

II 1/2 

into account, with in-plane and out-of-plane vortex spacings given by ajj = a^je^ 
and = e^^[;&)a^ respectively, we have 

4' = (35) 

.±2 _ ( 



so that for B || a6, 

l|2 '^0 /Q^7^ 

^0 = -TTTV (37) 

Thus, compared to the orientation || c, a vortex in the a6-plane can wander a similar 
distance in the plane, but a much smaller distance perpendicular to it. The trapping 
area i^\u^ = \e\j e\){^y\u^^., so that the length of a pinned vortex segment ^c" = \jnii}\u^ 
exceeds £ for B || c by a factor e\){^)je\. The critical currents for vortex motion within 
the ab plane and perpendicular to it are 

3f = ^n,ulu^ = -^f^ (B II ab;j \\ c; B > BD (39) 

r_\_ 2 

iat = i-ni4u^ = -^^Jt; (B II a6; j || ab.j ±B;B> B;) (40) 
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The crossover field B^= /si^^O) delimiting the range of vahdity of Eqs. ( |39[|40| ) with 
respect to Eqs. ( 27p8 ) is determined by comparing the high- and low field values of 
the respective critical current densities. It depends only on the specific details of the 
pinning centers. 

The ratio j = of the in-plane and out-of-plane critical current densities 
is once again determined by the coherence length anisotropy. Note that where, in the 
weak pinning regime, the critical current density for B || a6, j || a6, is the same as (or 
comparable to) that for B || c, j || a6, in the strong pinning regime it may be noticeably 
smaller. 

Finally, we give the expression for the field-dependent a6-plane critical current 
density as function of field orientation, 

valid for fields lower than the angle-dependent crossover field = B"" je\){p\ For 

small spherical defects in a single-band superconductor, this yields the well-known 
scaling law in e^B 




Figure 4. (Color online) Angular dependence of the critical current density in 
the single vortex limit of strong pinning, for anisotropics representative of (a,b) the 
(Ba,K)Fe2As2 system and (c,d) the NdFeAsO and SmFeAsO systems, for ellipsoidal 
oblate and prolate insulating defects with half axes = 5^, 6^ = ^, and bx = 
bz = 5^, respectively, is the angle with respect to the a^plane. 
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3. Discussion 

The following salient features emerge from the above derivations. First is the field- 



angular dependence of the critical current density in the a6-plane. Eqs. (26) and (41) 
for the strong pinning critical current both involve the ratio of powers of the penetration 
depth anisotropy and the coherence length anisotropy. The penetration depth anisotropy 
in the denominator mainly comes from the line energy term appearing in the elementary 
pinning force. The coherence length anisotropy appears mainly due to the varying 
core size for different orientations, which determines the order parameter gradient and 
therefore, also, the magnitude of the pinning force. 

As a result, for small defects of size comparable to the coherence length, only the 
latter contribution matters. Fig. [2](a,c) shows the field-angular dependence of the ab- 
plane critical current density for such small defects, for two values of the penetration 
depth anisotropy, and different e^. Such a plot mimics the temperature dependence of 
pinning in the iron-based superconductors, in which < Sx increases with temperature, 
the two reaching near equality near Tc. At low temperature <C and the coherence 
length anisotropy completely dominates the behaviour. This is why thin films and 
crystals of the 1111-family of iron-based superconductors seem to behave as single-band 
anisotropic superconductors, and scaling of jah with the product of the magnetic field 
and Eq. ([T]) yields the Sc2-anisotropy ^23^. Note that the change of apparent anisotropy 
with changing may be very difficult to observe; for example, as the temperature is 
increased, the absolute value of the coherence length also increases, so that defects that 
are to be considered "large" at low temperature become "small" close to T^. 

For larger defects, see Fig. [2](b,d), the scaling is disrupted, since now, the role 
of the angular-dependent vortex line energy becomes important. This results in a 
local maximum of the critical current for field aligned with the c-axis. For sufficiently 
large or strongly-pinning defects, this effect can result in the inversion of the apparent 
material anisotropy. As summarized in Fig. [3| this effect is not peculiar to multi-band 
superconductors, but is a property of strong fiux pinning, as opposed to weak pinning. 

The effect of large defects on the field-angular dependence of the critical current 
density is particularly pronounced in the case of oblate and prolate defects, see Fig. [4} 
Note that this effect is different in nature from that of extended columnar [34l [35l [36] , 
planar defects [l5l US] or intrinsic pinning [HI HHl HHl ED], for which depinning takes 
place by half- [35l [47] or quarter-loop [5l] nucleation rather than by pin-breaking. The 
nucleation-type depinning leads to a sharp maximum of the critical current density 
for field aligned with the defect direction, while the depinning from oriented ellipsoidal 
defects leads to a much broader maximum. In Fig. [4j the broad maxima of the critical 
current density as function of field angle correspond to the effect of the ellipsoidal defects, 
while the sharp structure, for field || a6 is due to the material anisotropy. On top of that 
comes the possible effect of intrinsic pinning by the layered structure of the material, 
such as this may have been revealed by deviations from field angular-dependent scaling 
in Ref. [23j. The above analysis shows that in presence of strong pinning, the breakdown 
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of simple scaling with can arise because of several reasons, among which, the 

presence of point-like defects of dimensions larger than the coherence length, and the 
presence of defects of distinctly different dimensions and aspect ratios. 

A final and general conclusion is that the measurement of the critical current density 
anisotropy for field parallel to the a6-plane (perpendicular to the anisotropy axis of the 
material) is a useful tool for the direct determination of the coherence length anisotropy 
It was shown above, on very general grounds, for both the weak- and strong pinning 
scenarios, that the jf^/j^ ratio is always equal to e^. This is because the relative 
pinned volume, i.e. the collective pinning length in the case of single vortex collective 
pinning, the correlation volume Vc in the case of bundle pinning, and the pinned length 
C in the case of strong pinning, only depends on the field orientation and not on the 
direction in which the vortices are being driven. Moreover, the critical current density 
does not depend on the range of the pinning potential. Therefore, only the anisotropy 
of the pinning force determines, in the end, the anisotropy of the critical current density. 
Since this is given by the anisotropy of the vortex core in a given field orientation, we 
may expect that it directly yields the coherence length anisotropy. 

The latter fact was recently exemplified by measurement of the three critical 
currents j^^, j^^, and j^^ on single crystalline LiFeAs [26]. Figure [5] shows the 
temperature dependence of the three critical currents, as determined by Hall array 
magnetometry of the fiux density gradient on different crystal faces under different 
modes of field application. From this, one may determine, directly, the ratio jf^/j^ to 
yield the temperature- and field-dependent coherence length anisotropy factor 6:^(T, B). 
Fig.[6]shows that the result is in excellent agreement with the same quantity, determined 
from the upper critical field measurements of Refs.p!0l [Tn[52] . Note that the anisotropy 
factor has a non-trivial behaviour, in which it increases as function of field, but 
decreases as function of temperature. It can therefore not be simply explained in terms 
of the reduced field B / 3^2 at which the measurements were performed, in conjunction 
with the temperature dependence Bc2{T). 

4. Conclusion 

The field-angular dependence and anisotropy of the weak and strong pinning critical 
current density in uniaxial anisotropic single-band and multi-band superconductors was 
considered, using a simple phenomenological approach introducing distinct anisotropy 
factors Sx and for the penetration depth and the coherence length respectively. 
Expressions for the field-angular dependent and anisotropic critical current density 
are derived. It turns out that the field-angular dependence of the a6-plane critical 
current density is determined by the competition between two effects. These are the 
angular variation of the vortex line energy, and therefore the pinning energy, which 
is maximum for field parallel to the c-(anisotropy) axis, and the angular variation of 
the vortex core size (coherence length), which yields larger elementary pinning forces 
for field parallel to ah. As a result, the field- angular dependence of the critical current 
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r(K) 

Figure 5. (Color online) Temperature dependence of the three independent critical 
current densities along the main crystalline axes of single crystalline LiFeAs, at an 
applied field of 1 T. The data depict j^^ (o, current density along c, field along a6), j^^ 
(o , current in the a6-plane, perpendicular to the magnetic field, also aligned in the 
a6-plane), and j^^ ( A, current in the a6-plane, field along the c-axis). The relevant 
magnitude of j^^ and j^^ suggest pinning either by insulating point defects of typical 
size ~ ^ — 2^, or by areas larger than 2^ showing weaker superconductivity than the 
surrounding matrix. 

density in materials with a relatively large coherence length anisotropy (small e^) and/or 
relatively small point-like pinning centers is indistinguishable from that in a single-band 
superconductor. Analysis of this angular dependence then yields the coherence length 
(upper critical field) anisotropy. On the other hand, in (single-band and multi-band) 
superconductors with modest anisotropy and large point-like pinning centers, strong 
pinning for field oriented along the c-axis can lead to an apparent inversion of the 
anisotropy. Furthermore, we have shown, on very general grounds, that the ratio of 
the c-axis and a6-plane critical current densities for field along the a6-plane always 
yields the coherence length anisotropy e^. This can therefore be extracted as function 
of temperature and field, at magnetic inductions far below the upper critical field Bc2- 
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